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SUMMARY 


Improving the aerodynamic characteristics of an 
airplane with respect to maximizing lift and minimizing 
induced and parasite drag are of primary importance in 
designing lighter, faster, and more efficient aircraft. 
Previous research done by Olson (Ref. 18) has shown that a 

properly designed biplane wing system can perform superiorly 
to an equivalent monoplane system with regard to maximizing 
the lift to drag ratio and efficiency factor. Biplanes 
offer several advantages over equivalent monoplanes, such as 
a 60% reduction in weight, greater structural integrity, and 
increased roll response. The purpose of this research is to 
examine, both theoretically and experimentally, the 
possibility of further improving the aerodynamic 
characteristics of the biplane configuration by adding 
winglets. Theoretical predictions were carried out 
utilizing vortex-lattice theory, which is a numerical method 
based on potential flow theory. Experimental data were 
obtained by testing a model in the Pennsylvania State 
University's subsonic wind-tunnel at a Reynolds number of 
510,000. Results indicate that the theoretical predictions 


agree fairly 

well 

with 

the experimental results. 

More 

importantly , 

the 

results showed 

that the 

addition 

of 

winglets improved 

the 

performance 

of the 

biplane 

with 


respect to increasing the lift-curve slope, increasing the 
maximum lift coefficient, increasing the efficiency factor, 
and decreasing the induced drag. 
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coordinates of a particular horseshoe vortex 
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CHAPTER I 


INTRODUCTION 


In modern times, with the advent of wide-body jets and 
supersonic transports, the aircraft has proven to be a very 
fast and efficient means of transportation from a 
pa s senge r-sea t-mi le s-pe r-ga 1 Ion standpoint. However, 
private, corporate, and business travel in light aircraft, 
commonly referred to as general aviation, also plays a vital 
role in America's transportation system. For transportation 
needs between cities located approximately 100 to 500 miles 
apart, the light aircraft is a very attractive means of 
transportation. For example, a small four-place 
single-engine aircraft, such as a Mooney 201, will transport 
four people at 187 miles per hour while burning 9.9 gallons 
of fuel per hour. This translates to 18.9 statute miles 
traveled per gallon of fuel burned, or 75.6 
passe nge r-s ea t-mi 1 e s-pe r-ga 1 Ion . This is not as efficient 
as a small automobile. However, when considering the high 
speed and straight line travel which the airplane affords, 
it appears to be quite advantageous. Compare this to a 
McDonell-Douglas DC-10 wide-body jet, which can transport 
255 passengers at a specific range of 0.22 miles traveled 
per gallon of fuel burned. This yields 56.1 
passenger-seat-miles-per-gallon. From this standpoint, it 
can be seen that the general aviation airplane is the more 
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efficient means of transportation. However, for the general 
aviation airplane to remain an efficient means of 
transportation, continous . improvements must be made on 
improving aerodynamic efficiency and reducing drag. The 
high cost of aviation fuel, coupled with increasing 
operating and maintenance costs, is forcing aircraft 
manufacturers to design and build more efficient and better 
performing aircraft. The National Aeronautics and Space 
Administration (NASA) has been performing a great deal of 
research in the area of improving aerodynamic efficiencies. 
Such research projects include experimentation with winglets 
and wing-tip extensions on general aviation aircraft, as 
well as first and second generation jet transports. It has 
been determined that properly designed winglets can 
significantly reduce induced drag at cruise lift 
coefficients without imposing severe additional structural 
loads ( Re f • 4 ) . 

It is the purpose of this research to study, both 
experimentally and theoretically, the aerodynamic 
characteristics and aerodynamic efficiencies of a biplane 
configuration utilizing winglets. In present literature, 
there is very little material available on the aerodynamic 
theory of biplanes. In the early days of aviation, when the 
monoplane was first introduced, research on biplane theory 
was virtually discontinued. The analytical determination of 
biplane characteristics was very complicated due to the 
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complex interactions between the two wings operating in 
close proximity to each other. Until recently, very few 
attempts have been made to optimize the aerodynamic 
efficiency of the biplane. 

If a biplane wing system could be designed to operate 
as efficiently as an equivalent monoplane system (with the 
same equivalent wing loading and aspect ratio), the biplane 
would offer several advantages. Because of the decreased 
structural constraints of the biplane, the biplane wing 
system can be as much as 60% lighter in weight than the 
equivalent monoplane system. Also, because of the increased 
roll response of the biplane, much less aileron area is 
required. This means that most of the wing’s trailing edge 
can be utilized for high lift devices, such as fowler flaps • 
Biplanes offer the potential for excellant low-speed 
maneuverability, good short-field performance, good load 
carrying capability and rugged construction. 

Previous Investigations 

There are three terms commonly used to define the 
geometry of a given biplane configuration. They are gap, 
stagger, and decalage. The gap(Ga) is the distance one wing 
is located above the other measured in percent chord length. 
Stagger(St) is the distance the upper wing is ahead of, 
positive, or behind, negative, the lower wing measured in 
percent chord length. De ca 1 ag e ( De c ) is the angle between 
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the chord lines of the upper and lower wings. The decalage 
angle is negative V7hen the lower wing is at a greater angle 
of attack than the upper wing. 

In 1918 F.H. Norton (Ref. 17) conducted experiments 
utilizing three-dimensional non-symme tr ical biplane 
airfoils. His results showed that maximum aerodynamic 
effiency is achieved at the highest degree of stagger 
physically possible. He varied only the stagger while 
holding the decalage constant at 0 degrees and the gap 
constant at one chord length. Also, Norton discovered that 
positive stagger greatly reduces the center of pressure 
travel, which simplifies the problem of stability. 

In 1929 Knight and Noyes (Refs. 11-13) conducted 
several three-dimensional non-symmet r ical biplane airfoil 
tests and concluded that increasing stagger in the positive 
direction, or increasing the gap, tends to equalize the 
loads on the two wings (this does not entirely agree with 
the results predicted theoretically by the vortex-lattice 
computer program). They also discovered that changes in 
decalage from 0 degrees for the orthogonal biplane (stagger 
equal to zero and gap equal to one chord length) tended to 
reduce the maximum lift coefficient. This is because (for 
the orthogonal case) the greatest maximum lift coefficient 
is reached when both wings stall nearly together. This 
occurs when they are at the same effective angle of attack. 
If positive stagger is present, the lower wing must be 
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operating at a higher angle of attack (negative decalage 
angle) than the upper wing in order to have sufficient stall 
natch. This effect is due to the lower wing being emersed 
in the downwash of the upper wing, thereby reducing the 
effective angle of attack of the lower wing. 

In 1936 r[ . Nenadovitch (Ref. 16) conducted several 
experiments to determine the aerodynamic characteristics of 
two-dimensional symmetrical biplane airfoils. He discovered 
that at a gap of one chord length, a stagger of one chord 
length, and a decalage angle of -6 degrees there was a 
substantial reduction in drag. 

These results are significant. However, in none of the 
Previous experimentation has any comparison been made to an 
equivalent monoplane configuration until 1974. In 1974 E.C. 
Olson (Ref. 18) conducted extensive experimentation on 
three-dimensional non-symmetrical airfoil biplane 
configurations in which the geometry was varied about 
Nenadovi tch ' s optimum test configurations. However, in 
conjunction with the various biplane configuration tests, 
Olson also tested an equivalent monoplane system. This is a 
monoplane system which has the same wing area as the biplane 
system as well as a similarily related aspect ratio. He 
discovered that at certain optimized geometric 
configurations, the biplane outperformed the monoplane 
configuration with respect to minimum drag and maximum lift 
to drag ratio. The biplane configurations were tested with 
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and without the fuselage. Specifically, Olson's experiments 
resulted in the following conclusions: 


1. At a gap of one chord length, a stagger of 0.875, 

and a decalage angle of -6 degrees, the biplane 
configuration showed a 25% reduction in drag over 
the monoplane at a typical cruise lift 

coef f icient . 

2. At a gap of one chord length, a stagger of 0.875, 
and a decalage angle of -5 degrees, the biplane 
configuration showed a 31.2% increase in the 
maximum lift to drag ratio while producing a 21.4% 
reduction in drag over the monoplane. 

3. The most effective overall biplane configuration 
was found at a gap of 0.875, a stagger of one 
chor<^ length, and a decalage angle of -6 degrees. 
This biplane configuration showed a 16.3% increase 
in the in the maximum lift to drag ratio, and a 
14.3% reduction in drag at a lift coefficient of 
0. 175. 

4. All biplane configurations showed a reduction in 
the maximum lift coefficient when compared to an 
equivalent monoplane configuration. 


Past research has shown that the addition of winglets 
to a given wing configuration can significantly reduce the 
induced drag. The presence of the winglets causes a 
physical constraint to the flow field near the location of 
the winglet, which is usually at or near the wing-tip. This 
constraint weakens the strength of the trailing vortices 
shed near the vicinity of the wing-tip. This reduction in 
the strength of the trailing vortices causes a reduction in 
the induced downwash, particularly in the vicinity of the 
outboard section of the wing. By reducing the induced 
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down wash, the effective angle of attack of the wing is 
increased. This results in a more even spanwise load 
distribution across the wing; the net result being a more 
efficient wing. 

The geometric configuration of the winglet is primarily 
described by two parameters: the winglet cant angle and the 
winglet toe angle. The cant angle is defined as the angular 
deflection of the winglet planform relative to a vertical 
plane which is perpendicular to the aircraft's lateral axis. 
At a cant angle of 90 degrees, the winglet acts as a 
wing-tip extension, which is unfavorable because of the 
increased bending stresses imposed on the wing structure. 
Also, a cant angle of 0 degrees is unfavorable due to the 
increase in interference drag caused by thickening boundary 
layer interactions at the wing-winglet joint. The toe angle 
is the incident angle of attack at which the winglet is 
mounted relative to the airplane's longitudal axis. 

In addition to altering the spanwise. load distribution, 
the winglet can also induce a negative drag contribution. 
This is caused by the forward tilting of the winglet normal 
force vector. This forward tilting effect is caused by the 
winglet operating at an induced angle of attack, which is 
brought about by the vectorial addition of the sidewash 
velocity and freestream velocity vectors. 
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Description of Research 

The purpose of this research is to investigate the 
possibility of further increasing the aerodynamic 
performance of the biplane configuration by adding winglets 
to the already optimized biplane configurations found in 
Olson's experiments. Pertinent aerodynamic characteristics 
of the bi pi ane-wi ngl e t configuration will be predicted 
theoretically as well as determined experimentally. If the 
biplane- winglet configuration could be optimized to the 
extent that it could perform as well as an equivalent 
monoplane (with respect to minimizing drag and maximizing 
the lift to drag ratio), the biplane could offer several 
advantages, which have been previously mentioned. 

The first step in this research was to theoretically 
analyze and predict the aerodynamic performance of a given 
biplane configuration with and without winglets. The method 
used was a finite-element, three-dimensional potential flow 
code , commonly refered to as the vortex-lattice method. 
Vortex-lattice utilization is commonly used throughout 
industry and government research to predict subsonic 
aerodynamic characteristics of complex planforms as well as 
predicting spanwise and chordwise load distributions on 
aerodynamic structures. Research has shown vortex-lattice 
theory to predict aerodynamic characteristics of complex 
planforms with considerable accuracy. 
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The second step was to experimentally test an already 
optimized design configuration. Based on previous research 
done by Nenadovitch and Olson, it was determined that the 
optimum test case would be a biplane-wingle t configuration 
with a gap of one chord length, and a stagger of one chord 
length. The configuration was tested with and without 
winglets at a Reynolds number of approximately 510,000. The 
decalage angle was varied from 0 to -5 degrees. The 
experimental data was then reduced and several wind tunnel 
correction factors applied to yield experimentally correct 
lift and drag data. 
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CHAPTER II 

PRESENTATION OF THEORY 

Before the advent of modern high-speed computers, 
simple problems in aerodynamics had to be solved using 
classical theory. As the problems became more complex, the 
application of classical theory became quite cumbersome due 
to the complex conformal transformations which had to be 
utilized. It was very difficult to optimize a given design 
except through trial and error. The analytical approach to 
predicting the aerodynamic characteristics of a relatively 
simple biplane configuration proved to be quite cumbersome, 
even for a simple mathmatical model such as the classical 
bound-vortex lifting-line method. The addition of winglets 
to the biplane configuration makes the problem much more 
complex. The interference effects of all components must be 
considered, since they have a significant effect on the 
induced drag and spanwise load distribution. 

Basically, * the vortex-lattice method is a 
f ini t e- e lement method which utilizes a vortex-lattice 
representation of the aircraft's lifting surfaces coupled 
with classic equations and theorems for computing 
aerodynamic characteristics such as lift, induced drag, 
spanwise load distributions, and wing efficiency factors. 
This method assumes steady, irrotational , inviscid, 
incompressible, attached flow. Therefore, numerical results 
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can only be assumed valid at subsonic speeds when the wing 
system is operating at a less than critical angle of attack. 
The vortex-lattice method is commonly used for predicting 
the aerodynamic characteristics of complex three-dimensional 
planforms such as the Lockheed boxplane and Whitcomb winglet 
configuration (Ref. 4). In this research, the 
vortex-lattice method will be utilized in predicting the 
aerodynamic performance of the biplane-winglet 
configuration. 


Basic Theoretical Concepts 

Fundamental to the development of the vortex-lattice 
model is the representation of the aircraft's non-planer 
lifting surfaces by a system of rectangular horseshoe 
vortices. Basically, each planform surface is divided into 
several finite elemental panels which extend chordwlse and 
spanwise across the entire planform surface. At the quarter 
chord point of each elemental panel a bound horseshoe vortex 
is located, and at the three-quarter chord point a 
corresponding control point is located. Figure 2-1 shows a 
typical section of wing which has been broken down into 
several elemental panels. At each control point the no-flow 
through condition must be satisfied; that is, the flow must 
be tangential to the planform surface at this point. This 
concept first appeared in a paper by E. Pistolesi 
1937 (Ref. 5). He found that by using the 1/4-3/4 chord 
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rule ~ section lift and moment predictions for a cambered 
airfoil at a constant angle of attack were exactly that of 
thin-airfoil theory. In 1942 J. Weissinger applied this 
method to wing configurations of finite aspect ratio and 
also achieved accurate results. This method has been widely 
accepted and is used throughout present research which 
utilizes vortex-lattice methods. 

The first step in utilizing v o r t ex- la t t i c e theory, 
assuming the planform geometry has been defined, is to 
determine the number of chordwise and spanwise horseshoe 
vortices that are to be located on the planform surfaces. 
The three-dimensional coordinates which locate the bound 
horseshoe vortices are next computed, as well as computing 
the coordinates of the coo re s pond ing control points. Once 
this is accomplished, the induced velocities from the total 
vortex system can be equated to the freestream velocity 
component normal to the lifting surface at each control 
point. Application of the tangent flow boundary condition, 
(assuming a symmetrical loading) will yield a set of N 
siir.ul tanious equations, each equation consisting of N 
unknown horseshoe vortex strengths. The fundamental laws of 
induced velocity from a vortex filament are utilized in 
calculating. the horseshoe vortex induced flow-field at each 
control point. Once the set of N simultanious equations is 
solved, and the strength of each individual horseshoe vortex 
is determined, the Ku t t a- J o uko ws ki theorem for lift from a 
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vortex filament is utilized to determine the section lift 
coefficient. Finally, the finite wing lift coefficient can 
be obtained by numerically integrating the spanwise load 
distribution across the entire planform surface. 

The induced drag created by the bound vortices located 
on the planform surfaces, which is of primary importance, 
can be determined for any given loading and operating 
condition by utilizing the following basic laws and 
theorems: Biot-Savart Law, Kut ta- Joukowski Theorem, and 
flunk's Theorems I and II. 

flunk's first theorem (Ref. A) can be stated as follows: 

The total induced drag of any multi-plane lifting 
system is unaltered if any of the lifting elements 
are moved in the direction of the motion provided 
that the attitude of the elements is adjusted to 
maintain the same distribution of lift among them. 

This theorem is commonly referred to as "Munk's stagger 
theorem." An illustration of this theorem is shown in 
figure 2-2. Several practical applications can be reasoned 
from this theorem. First, the chordwise distribution of 
pressure does not affect the theoretical induced drag of the 
aircraft if constant section lift is maintained. Second, 
wing sweep or biplane stagger does not affect the 
theoretical induced drag as long as the spanwise 
distribution of lift is constant. A third application is 
that the load from a system of multi-surfaces, such as the 
wing and horizonal tail, with the same projection in the Y-Z 
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plane can be made equivalent to a single surface for the 
purpose of calculating induced drag. 


In the 

following theoretical 

development , 

use 

will 

be 

made of flunk's first theorem 

to combine 

the 

chordwise 

distribution 

of vorticity into a 

single chordwise 

load 

and 

to translate 

all loads into the 

0,Y,Z plane. 





flunk's second theorem is illustrated in figure 2-3 and 
can be stated as follows: 


In calculating the: total induced drag of a lifting 
system, once all the forces have been concentrated 
into the 0,Y,Z plane, one may, instead of using 
the actual values of the velocity normal to the 
lifting elements [ V Cx»y,z) ] at the original 
points of application of the forces, use one-half 
of the limiting value of the normal velocity 
[ 00 '* y * z ) 1 f° r t * le corresponding values at 
points P(0,y,z). 


This theorem allows the computations to be done in the 
Trefftz plane, a plane which is located infinitely far 
downstream, rather than in the real plane. In the 
subsequent theoretical derivation, this fact will be 
utilized in order to make all the induced drag computations 
in the Trefftz plane, thereby greatly simplifying the 
calculations. 


The third theorem given by Munk is presented as 
follows : 


When all the elements of a lifting system have 
been translated longitudinally to a single plane, 
the induced drag will be a minimum when the 
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component of the induced velocity normal to the 
lifting element at each point is proportional to 
the cosine of the angle of inclination of the 
lifting element at that point. 

This theorem is illustrated in figure 2-4 and can be 
summarized in equation form as: 

V = w cos 0 2.1 

n o 

For a horizonal lifting element it can be seen from equation 
(2.1) that the normal velocity (downwash) across the span is 
equal to a constant. For a vertical plane ( 6 =90 degrees), 
the normal velocity (sidewash) must be equal to zero for 
minimum induced drag. The physical interpretation of this 
theorem will be further illustrated in a subsequent section. 

The basic equation for calculating the induced drag can 
be derived by applying the Kut ta-Joukowsaki theorem in the 
drag direction. By utilizing Munk's theorems, the 

calculations can be accomplished in the Trefftz plane. 
Thus, the basic equation for calculating the induced drag 
for an arbitrary non-planer lifting system, expressed in 
terms of the Trefftz plane and using vector notation is: 

D i * 2 Y - / V • n N di, 2.2 

00 

This integral is a line integral taken around the 

perimeter of the projection of the lifting surface in the 
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Trefftz plane. The vector V is the resultant induced 
velocity vector in the Trefftz plane from all horseshoe 
vortices located on the load perimeter. N represents the 
load perimeter normal force per unit span. For a horizonal 
lifting surface, K would represent the section lift force. 
The vector n is a unit vector, normal to the load 
perimeter. 

Physical Interpretation of Theoretical Concepts 

To provide a better understanding of induced drag 
calculations, the theoretical concepts discussed in the 
previous section will be illustrated using a monoplane 
wing-winglet configuration. In figure 2-6 the sources of 
induced drag for a wing-winglet configuration are shown. 
They are : 

Induced drag due to the induced flow by the wings 
on the wing 

Induced drag due to the induced flow by the wings 
on the winglet 

Induced drag due to the induced flow by the 

winglets on the wing 

Induced drag due to the induced flow by the 

winglets on the winglet 

For simplicity, the effects of symmetry are included in 
the sources of induced drag shown and are not delineated 
separately. 



In figure 2-5a the effect of the wing induced flow is 
shown. The wing under positive load produces a downwash on 
itself which results in the wing force vector, F , tilting 
rearward by an angle . The w'ing force vector is 
perpendicular to the resultant velocity vector, V , by 
definition of the Ku t t a- J oukows ki theorem. A sidewash is 
also produced by the wing at the winglet. As can be seen in 
figure 2-5b, the sidewash at the winglet combined with the 
freestream velocity produces a tilt forward of the resultant 
winglet force vector. This produces a negative drag 
component as well as a side force component on the winglet. 

In figure 2-5c, the induced drag resulting from the 
sidewash induced by the winglet on itself is presented. 
This results in a rearward tilting of the resultant force 
vector, which creates an attendant induced drag component as 
well as an additional side force component. It should be 
noted that the direction of the winglet force vector is 
consistant with a positive (upload) on the wing. The 
winglet also induces an upwash on the wing. In figure 2-5d 
it can be seen that this upwash rotates the resultant wing 
force vector forward. This produces an additional lift 
force on the wing as well as a negative drag component. 

The results of figures 2-5a, 2-5b, 2-5c, and 2-5d are 
summarized in figure 2-6, where all the induced velocities 
are combined. For minimum induced drag, equation (2.1) 
indicates that the velocity normal to the winglet must be 
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equal to zero. This can be seen to occur when the sidewash 

’ * .. . i 

produced on the winglet by the wing exactly cancels the 
sidewash produced by the winglet on itself. In other words, 
the induced angle of attack of the winglet is zero. The 
induced drag of the wing is also minimized by the presence 
of a winglet since the winglet causes a reduction in the net 
dowriwash at the wing; hence, the induced angle of attack is 
reduced. Also, the winglets allow the wing to be loaded 
more heavily out towards the tips, which of course results 
in a more efficient wing. 

Theoretical Application 
The Vo r t e x- La t t i c e Computer Program 

The purpose of this section is to present and discuss 
the methods used in developing the vortex-lattice computer 
program. Relevant equations and formulas will be discussed 
in order that the reader may understand, and if necessary, 
modify the existing program. Figure 2-7 shows the computer 
program flowchart which represents the internal structuring 
of the vo r t ex- 1 a 1 1 i c e computer program. 

The first function performed by the program is the 
input of the bi pi a ne-wi ng 1 e t planform geometry, as well as 
the input of data representing a given flight condition. 
Upon execution of the program, the computer will prompt the 
user for the following: biplane aspect ratio, biplane 
sta 8ger, decalage angle, wing twist, maximum winglet toe 
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angle, and angle of attack. The wings and winglets are 
assumed to have no camber or taper, and the biplane gap is 
held constant at one chord length. Since the flow is 
inviscid and incompressible, the velocity may be arbitrarily 
chosen. 

The next function performed, after the total number of 
spanwise and chordwise horseshoe vortices have been 
determined, is the computation of the coordinates of all the 
horseshoe vortices and their cooresponding control points. 
In this program, the number of chordwise and spanwise 
horseshoe vortices is pre-set. The upper and lower wings 
each contain a total of 80 horseshoe vortices (40 located 
spanwise and 2 located chordwise for each wing), while the 
winglets each contain a total of 16 horseshoe vortices (8 
located spanwise and 2 located chordwise for each winglet). 
These numbers were arrived at by analyzing output data from 
several program executions , . and choosing numbers which would 
yield converging results without consuming enormous 
computational time. 

In constructing the vortex-lattice planform model, only 
a half-span model is constructed since a symmetrical loading 


is assumed 

. The program computes the 

posi tion 

o f 

each 

horse shoe 

vortex, represented by pn^ 

, qn , m, 
n n 

where 

the 

subscript 

th 

n represents the n 

horseshoe 

vortex. 

Similar ily 

, the program also computes 

the position of 

each 


corresponding control point, represented by P v v > < l v v » rv v » 
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where the subscript v represents the position of the v^. 
control point. This is accomplished by applying the 1/4-3/4 
rule to each elemental panel. Figure 2-8 shows the 
half-span planform model of the b i p 1 a n e -wi ng 1 e t 
configuration represented by a system of rectangular 
horseshoe vortices. In this model, through ate 
located on the upper wing, through Tgg are located on 
the lower wing, and r Q , through r o , are located on the 
starboard winglet. The P,Q,R axis system is also shown 
with the origin located at the mid-span point of the leading 
edge of the upper wing. 


Since a symmetrical loading is assumed, the tangential 
flow boundary condition will be applied only to the control 
points located on the starboard planform. However, the 
total induced velocity at each control is contributed to by 
each and every bound horseshoe vortex located on both the 
starboard and port planforms. Because of the symmetrical 
loading assumption, T n on the starboard wing is equal to 
T n on the port wing. Therefore, only the half-span 
horseshoe vortex strength distribution must be solved for. 


Next, the coordinates, p , of the control point 

vn 

relative to the xi^ 1 horseshoe vortex can be computed in 

the X , Y , Z axis system. For the starboard wing: 

x = pv - pn 
vn r v n 


P = (x , ys , z ] 
vn vn vn vn-' 


ys = qv - qn 
J vn n v n n 


z ■ rv - rn 
vn v n 
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and for the port wing: 


P 

vn 



^vn’ 



x = pv - pn 
vn v n 


yp = qv + qn 
7t vn v n n 


z = r v - m 
vn v n 


The influence coefficients, which relate the induced 
velocity at the v control point to the strength of the n** 1 
horseshoe vortex (which is inducing the flow at that control 
point) must be computed next. They are computed in order to 
determine the vortex induced velocity at each of the control 
points located on the starboard planform. 


The downwash influence coefficients, which represent 
the induced velocities caused by the bound horseshoe vortices 
located on the starboard planform, can be computed from the 
following expression. The angle represents the angle of 
the bound horseshoe vortex filament in the P-R plane. For 
the horizonal surface (wing), <f> =0 degrees, and for the 
vertical surface (winglet), 0 =90 degrees. 


Fws 


-x cos d> 
vn in 


vn 2 .2 

x +(Z cos d> - ys sin <J> ) 
vn vn n J vn n 


(ys +s cos d> )cos 6 +(Z +s sin 6 )sin<t> 
_ vn Y n n vn Y n Y n 


[x 2 +(ys +s cos (f) ) 2 +(Z +s sin 4> ) 2 ] 2 
vn vn n vn Y n 



Page 22 


(ys - s cos 4> )cos d> + (Z - s sin <J> )sin <f> 
vn n n vn n n 


lx 


(ys - s cos )* + (Z - s sin <f> ) 2 
w vn r n vn r n 


vn 


(ys vn - s cos <> n ) 

2 2 
(ys - s cos 6) + (Z 11#t - s sin * ) 

vn n vn n 

0 * 

x 

i vn 

1 2 1 Y TV 

lx + ( y S xin " 8 COS O + ( Z ^ ” 8 Sin O J 2 
vn vn n vn n 

l J 

(ys vn + 8 cos <|> n ) 

(ys vn + 8 cos ^ n )2 + (Z vn; + 8 sin *n )2 

x __ X vn 

>*» + (ys vn + * c °“ V 2 + (Z vn + 6 sln V 2lV2 j 


Similarily, the downwash influence coefficients which 

represent the port wing can be expressed identically to the 

above equation except yp is substituted in place of 

vn 

ys . The downwash at the v fc ^ control point, induced by the 
^ vn 

n** 1 horseshoe vortex located on the port and starboard 
planform, represented by w vn » can expressed as: 


w 


vn 


r 

4? (1Vs vn + F "I’vn ) 


2.3 
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In a similar manner, the sidewash velocity at the 
control point induced by the n^ horseshoe vortex can be 
computed. The sidewash influence coefficients representing 
the bound vortices on the starboard wing can be computed 
from the following ex pr e s s i on ( Ref . 15 ). 


Fvs 


-x sin 4> 
vn n 


vn x 2 + (Z cos <p - ys sin $ ) 2 
vn vn n vn n 


(ys +s cos <j> )cos <j> + (Z +s sin <(> )sin <p 
vn n n vn n n 


[X 2 +(ys +s cos <j) ) 2 +(Z +s sin <t> ) 2 ] 2 
vn J vn n vn n 


(ys^_ - s cos * )cos <f> +(z, - s sin <j> )sin <j> 

vn n n vn n n 


[x 2 +(ys -s cos ) 2 +(Z -s sin <j> ) 2 ] 2 
1 vn vn Y n vn T n 


(Z vn - s sin » B ) 


2 2 
(ys - s cos d> ) + (Z - s sin 4> ) 

J vn T n vn T n 


1 - 


vn 


[x vn +(ys vn’ 


s cos 0 )^+(Z -s sin <}i ) y i 2 
Y n vn n J 


Z + s sin <|> 
vn n 


(ys + s cos <fr ) 2 +(Z +s sin * ) 2 
vn n vn n 


i - 


vn 


t X L + ^y S u« +S COS +S S±n 

vn vn n vn n 
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The side wash influence coefficients representing the port 

planform can be expressed as above except yp^ is again 

substituted in place of ys^ '• Similar, ily, the side wash 

velocity at the v *"* 1 control point induced by the 

horseshoe vortex located on both the port and starboard 

planform, represented by v , can be expressed as: 

vn 


v = ~ (F vs + Fvp ) 

vn Ait vn vn 


2.4 


The total downwash at a given control point is equal to 
the sum of the induced downwash contributions from each 
horseshoe vortex located on the entire planform surface. 
The total downwash at the control point can be computed 

from the following: 

N N r-r 

w = 7 w = T 7” (Fws + Fwp ) 2.5 

v L . vn L . 4 tt v vn vn-' 

n=l n=l u -> 

Similarily, the total sidewash at a given control point can 
be expressed as 

n n p r 

v = y v ■ T fFvs + Fvp M 2.6 

V n=l vn n=l L 47r Vn vn J 


Next, the tangential flow boundary condition at each of 
the control points located on the starboard planform must be 
satisfied. First the horizonal (wing) planform surface will 
be considered. By equating the freestream velocity and 



local angle of attack to the local downwash velocity, the 


following boundary condition for the control point can 

be formulated. 


M 


l 

v=l 


a 

v 



N 

l 


n=l 


w 


vn 


2.7 


This can also be expressed as: 


M 

l 

v=l 


a = 


4ttV 


M 

I 

V=1 


N 

I 

n=l 


T (Fws + Fwp ) 
[n> vn r vn ; J 


2.8 


where represents the local angle of attack of the 
control point* M represents the number of horseshoe 
vortices and cooresponding control points located on the 
horizonal planform surface, and N represents the number of 
vortices and control points located on both the horizonal 
and vertical planform surfaces. Expanding equation ( 2.8 ) 
yields M linear equations; each equation containing N 
unknown circulation strengths. They can be expanded to the 
f ol lowi ng : 


r i (fWs l,l + F “ P l,l )+ r 2 (Pus l,2 + n *l,2 )+ 
r i (Fws 2,l +rWl, 2,l )+ r 2 (Fus 2,2 + PwI, 2,2 )+ 


+ r N (Fw6 l,K +Fwp l,K>* 4 ’Vl 


F 1 (F ” S M, 1 + ^M, 1 )+ r 2 ' tFBs M,2 + F " I> M,2 )+ ' ' • + 


I V F ” S M,N +F '' I 'm,n' 


4itV a 

oo M 
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The remaining linear equations can be derived by 
applying the same tangential flow boundary condition to the 
vertical surfaces. This boundary condition can be written 
as: 


N 

l 

v=M+l 


_1 

V 


N 


l 

v=M+l 


N 


l 

n=l 


v 

vn 


2,9. 


where i v represents the incident angle of attack (toe angle) 
of the winglet. This can also be written as: 


N 

l 


V5=M+1 


i 

v 


1 

4ttV 

CD 


N N 

l l 


v=M+l n=l 



+ Fvp ) 
r vn-'_ 


2.10 


Expanding this yields the remaining linear equations 
necessary to determine the unknown circulation strengths. 

r i (FvS Mfl,l + FvP Mfl,l )+r 2 (Fvs Mfl,2 + Fvp M+l,2^ + ' * * + r N (FvS tffl,N + Pv:p MH,H* 


4ttV 


■Sh-i 


r i ^ FvS M+ 2 , 1 + Fvp M+2 , 1 )+ r 2 (FvS Mf 2 , 2 + Fvp Mf 2 , 2 )+ * * * + F N (FvS tff 2 ,N + FvP Mf2 ,1^ 


4ttV 


1 Mf2 


F 1 (FVS N . 1 + FvP N . 1 )+ F 2 (FVS N , 2 + , 2 )+ * ‘ ‘ + r N (FVS N ,N + ,N ) “ 4irV -S 


Finally, the unknown circulation strengths can be 
computed by s imul taniously solving the complete set of 
linear algebraic equations. This can easily be done by 
first expressing the equations in matrix form as: 


[A]{r n } - {b} 
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where [A] is the coefficient matrix, and {B} is the boundary 
condition matrix* The solution can be obtained by inverting 
the coefficient matrix and multiplying it by the boundary 
condition matrix. 

i r n } - [a] _1 {b} 


The wing and winglet section loading coefficients can 
now be computed from the known circulation strengths using 
the following well-known relationship. 


2T 


c V 


00 


2.11 


The total lift produced by the upper and lower wings can be 
expressed as : 

b/2 

L - 2 / p T(y) dy 2.12 

o 

Expressing this in lift coefficient form yields: 



b/2 


V S 1 

oo v O 


r(y) dy 


2.13 


Converting this into a numerical integration form results 
in : 



V S . 
00 w n=l 


T 2s 
n 


2.14 
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Similarily, the side force coefficient for the winglet can 
be expressed as: 


N 

l 


L VS. _ . 
vZ «> wZ n=Mfl 


r 2s 

n 


2.15 


The induced drag coefficients are computed next b-y 
first calculating the total vortex induced velocities, 
normal to the control points on the load perimeter in the 
Trefftz plane. This is done by first calculating the 
position coordinates of the v s control point (in the 

Trefftz plane) relative to the horseshoe vortex as 

follows : 

x ■ - 00 
vn 


ys ■ qv - qn 
J vn v n n 


yp - qv 4 qn 
/r vn v n n 


z = rv - rn 
vn v n 


The influence coefficients are then recalculated. By 
knowing the circulation strengths of each horseshoe vortex, 
the down wash velocity normal to each control point on the 
horizonal planform in the Trefftz plane can be computed as 
well as computing the sidewash velocity normal to each 
control point on the vertical planform surface. Thus, by 
virtue of Munk's theorems and the Kutta-Joukowsi theorem, 
the induced drag coefficients can be computed. 
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The induced drag coefficient, , which is caused by 

w,w 

the wing inducing a downwash on itself, can be computed from 
the following expression: 


M 


I » r s 


D w,w V 2 S n=5-l v w n 
00 w 


2.16 


where w represents the downwash induced at the v 


th 


w 


control point by the wing. The induced drag 


coefficient, , which is caused by the winglet inducing 

w£ ,w 

an upwash on the wing, can be computed from the following: 


C n - ~£— l v T s 2.17 

u w£,w v S n®v«l v w£, n 
00 w 

where w represents the upwash induced at the v C ^ control 
V w£ 1 

point by the winglet. 


Similarily, the induced drag coefficients, and 

w,wZ 

C ., which are caused by the wing and winglet inducing a 

WZ yWZ 

sidewash on the winglet, can be expressed as follows: 


D 2 

wZ ,w Z VS n=v«M+l 

oo w 


N 

I v r s 
w t n 


2.18 


N 


I u r 

w M T 


^VyvZ S n=v«Mfl V w n 

oo v? 


2.19 


The total induced drag is then equal to the sum of the four 
induced drag coefficients. It is of interest to note that 
the wing induces a negative drag component on the winglet. 
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and similarily, the vinglet induces a negative drag 

/ 

component on the wing. The wing eff-iciency factor can be 
calculated from the following expression: 



2.20 


The theoretical drag polar can be now be computed. 
However, first the minimum profile drag and the incremental 
change in profile drag due to lift must be added to the 
induced drag term to yield the total drag coefficient* 


CL « CL + A C + C 
D D . D D, 

min p i 


2.21 


Finally, the three-dimensional lift-curve slope can be 
determined once the lift coefficient has been determined for 
different angles of attack. 


As was expressed previously, it is common in biplane 
theory to compute the wing efficiency factor utilizing the 
biplane's equivalent monoplane aspect ratio. The equivalent 
monoplane aspect ratio can be calculated from the following 
expression: 


EMAR = 


2 2 
V (1 + V) 

(y 2 + 2 ov y + y 2 ) 


bj upper wing span 


lower wing area 
upper wing area 


o from Figure 186 (Ref. 20) 


lower wing span 
upper wing 6pan 




P a g e 3 2 


INDUCED DRAG | STAGGERED = INDUCED DRA 5 | UNSTAGG£RED 


SPAN LOAD UNSTAGGERED STAGGERED 



Figure 2-2 Illustration of Munk's First Theorem 


REAL PLANE TREFFTZ PLANE 



Figure 2-3 Illustration of Munk's Second Theorem 
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Figure 2-6 


Combined Sources of Induced Drag for a 
Wing-Winglet Configuration 




Figure 2-7 Vortex-Lattice Computer Program Flowchart 
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CHAPTER III 

EXPERIMENTAL APPARATUS AND PROCEDURE 


Before conclusions can be drawn on the aerodynamic 
characteristics and performance of the biplane-winglet 
configuration as predicted by the vortex-lattice computer 
program, experimental data must be collected,' correlated, 
and analyzed to determine the validity of the theoretical 
predictions. 


Description of Apparatus 

For the purpose of obtaining experimental data, a 
half-span wind-tunnel model of the biplane-winglet 
configuration was constructed by the author and tested in 
The Pennsylvania State University's atmospheric 
closed-return subsonic wind tunnel. A half-span model was 
tested, which was mounted vertically in the wind tunnel, 
extending upward through the floor of the test section. The 
model was designed and constructed to have a stagger of one 
chord length, a gap of one chord length, and a decalage 
angle that could be varied from 0 to -6 degrees. The design 
of the model's geometric configuration was based on the 
already optimized design configurations found by Nenadovitch 
and Olson in previous wind tunnel experimentation utilizing 
biplane configurations. 
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The . two finite wings utilized in the model were 
constructed from solid mohogany using a router-assembly tool 
designed specifically for this purpose. The airfoil used 
was a NACA 0012, and the wing's cross-section profile was 
sanded to within 0.015 inches of the exact NACA 0012 airfoil 
specifications. The two wings, which are identical, have a 
chord length of 7.875 inches, a half-span of 19.75 inches, 
and a maximum cross-sectional thickness of 0.96 inches. 
This yields an aspect ratio of five, since the configuration 
tested represents a half-span model. 

Spanning the tips of the two wings is a constant-chord 
winglet. The winglet was constructed from pine and was 
sanded to form a thin 3% thick symmetrical airfoil. The 
winglet was held in place by four screws, and could easily 
be removed for conducting tests with and without the winglet 
attached. Also, by shimming, the winglet toe angle could be 
varied by as much as plus or minus 2 degrees, although zero 
incidence was used in the experimental tests. Zero 
incidence was chosen because the vortex-lattice computer 
program predicted that no aerodynamic advantages could be 
gained by having any winglet incidence present. The 
presence of a high-degree of toe angle incidence caused a 
slight increase in the induced drag. This was due to the 
increase in the spanwise load distribution along the 
winglet. The theoretical cases that were examined assumed 
the toe angle to be maximum negative where the winglet 
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joined^the upper wing, and maximum positive where the 
winglet joined the lower wing. The winglet was also assumed 
to have a linear twist between the points of maximum 
incidence. 

The chord length of the winglet was a constant 7.875 
inches, and no cant angle was present. Figure 3-1 shows the 
model mounted in the four by five foot test section of the 
subsonic wind tunnel. 

The bases of the wings were mounted, using 1/A inch 
bolts and angle brackets, to a 20 inch diameter disk plate, 
which was fabricated from 3/A inch laminated plywood. To 
vary the angle of attack of the model, the base disk was 
rotated about its geometric center. The aft wing was 
mounted on a 10 inch diameter disk plate which could be 


rotated 

about 

the aft 

wing's 

quarter 

chord 

point. This 

allowed 

the 

decalage 

angle 

to be 

varied 

from 0 to -6 

degrees 

. The 

winglet 

was designed to 

be used 

for decalage 

angles 

ranging 

from 0 

to -6 degrees . 

To reduce 

interference 


drag, the wi ng-wi ng le t joint was filleted with putty during 
experimental tests. Figure 3-2 shows the biplane-winglet 
model before being mounted in the wind tunnel. 

Located just below the wind tunnel test section is a 
six channel pyramid type strain guage balance which is used 

i 

for recording lift and drag forces during wind tunnel 
testing. The 20 inch diameter base disk, which supports the 
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model, was mounted on the pyramid balance prior to 
calibration of the balance for wind tunnel testing. 

Experimental Procedure 

Before the actual wind tunnel testing could begin, the 
wind tunnel balance had to be carefully calibrated. This 
was done by applying known forces on the balance in the lift 
and side force directions. The side force channel actually 
indicates the lift force on the biplane configuration. This 
is due to the non-standard method used in mounting the model 
vertically upright in the test section instead of 
horizonally across. Frictionless pulleys, nylon string, a 
level, and known weights were used to apply exact known 
loads to the wind tunnel balance. Several calibration tests 
were conducted and the following errors in the lift and drag 
measurments were determined; approximately 5-6% errror in 
the drag force measurments, and approximately 1-2% error in 
the lift (side) force measurments. The exact amount of 
error depends on the magnitude and range of the applied 
forces. The wind tunnel calibration correction factors were 
later applied in reducing experimental data. 

The biplane-winglet configuration was tested in the 
wind tunnel at a velocity of approximately 149 feet per 
second. The ambient air temperture in the test section 
varied from 115 to 120 degrees Farenheit while the 
barometric pressure was constant at 29.01 inches of mercury. 


This resulted in a Reynolds number of approximately 510,000. 

In the first set of experimental tests, the biplane 
configuration was tested with and without winglets at "-a 
decalage angle of 0 degrees, and a stagger and gap both 
equal to one chord length. The angle of attack was varied 
from -2 to 21 degrees in 2 degree increments. The angle of 
attack of the biplane configuration is actually the angle of 
attack of the upper wing, and the decalage angle is the 
incident angle of attack of the lower wing relative to the 
upper wing. 

In the second set of experimental tests, the same 
biplane configuration was tested using the same procedure 
except the decalage angle was changed to -5 degrees. The 
angle of attack was varied from -4 to 14 degrees, also in 2 
degree increments. 

Reduction of Experimental Data 

The first step in reducing the experimental data was to 
apply the wind tunnel balance correction factors (previously 
determined from calibration tests) to the measured lift and 
drag forces in order to yield the actual lift and drag 
forces. Next, standard wind tunnel correction factors were 
applied to yield experimental results that would be 
equivalent to the results obtained if the model were tested 


in f ree air 
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The first correction, known as 'horizonal buoyancy', 
was found to be negligible for the configuration tested. 
Buoyancy, which results from a thickening boundary layer on 
the test section walls, causes a decreasing static pressure 
gradient in the test section. 

The next correction, known as 'solid blocking', is 
caused by the physical constraint of the flow field normal 
to the flow direction. This results in a local dynamnic 
pressure increase over the model which tends increase the 
measured lift and drag forces. Also, due to physical 
constraints, the flow field surrounding the model's wake is 
also constrained. This effect, which is known as 'wake 
blocking', results in an increase in the measured drag 
force. Solid blocking and wake blocking have the same 
effect as horizonal buoyancy, which is an increase in the 
dynamic pressure over the model. 

Finally, a correction must be made to the angle of 
attack. This is also due to the physical constraint of the 
test section walls, which alter the trailing vortex system 
behind the wing. This causes a reduction in the effective 
angle of attack, which is caused by a reduction in the 
induced do wnwash. 

Tables 3.1 and 3.2 represent the data collected during 
wind tunnel testing after the wind tunnel correction factors 
have been applied to the experimental lift and drag data. 
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The wind-tunnel correction for angle of attack 
downwash constraint) was found to be negligible. 


( induced 


t 
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Table 3.1 


Reduced Experimental Data 
Stagger= 1.0, Cap= 1.0, Decalage= 0.0 

With winglet No winglet 

a 

C L. 


C D 


C D 

-2 

-0.093 


0.018 

-0.096 

0.016 

0 

0.000 


0.017 

0.008 

0.015 

2 

0.073 


0.018 

0.061 

0.016 

4 

0.220 


0.024 

0.201 

0.021 

6 

0.354 


0.036 

0.330 

0.033 

8 

0.519 


0.061 

0.487 

0.057 

10 

0. 654 


0.089 

0.6 28 

0.086 

12.5 

0. 775 


0.115 

0.726 

0.114 

14.5 

0.791 


0.156 

0.758 

0.148 

16.5 

0.827 


0.230 

0.784 

0.228 

18.5 

0.850 


0.267 

0.819 

0.265 

21 

0.772 


0.356 

0.753 

0.368 
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Table 3.2 

Reduced Experimental Data 



Stagger* 

1.0, Gap — 1.0, 

Defcalage 

* -5 


With 

wi ng le t 

No 

winglet 

a 

C L 

C D 

C L 

C D 

-4 

-0.0G1 

0.023 

-0.091 

0.021 

-2 

0.037 

0.022 

0.021 

0.020 

0 

0.186 

0.027 

0.220 

0.024 

2 

0.334 

0.038 

0.301 

0.035 

4 

0.4 68 

0.057 

0.425 

0.053 

6 

0.577 

0.075 

0.539 

0.072 

8 

0.698 

0. 104 

0.653 

0.100 

10 

0.791 

0.131 

0.744 

0.128 

12 

0.867 

0. 163 

0.833 

0.165 

14 

0.775 

0.260 

0.752 

0.256 



* I 


( « 



Figure 3-1 Biplane-Winglet Model Mounted in Wind-Tunnel Test Section 


Page 47 





Figure 3-2 Biplane-Winglet Model and Base Support 
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CHAPTER IV 

ANALYSIS OF RESULTS AND DISCUSSION 

Overall, the experimental results agreed fairly well 
with the theoretical vortex-lattice computer program, 
however, the experimentally determined drag coefficients 
were found to be quite large. This was discovered to be an 
effect caused by testing the biplane-winglet configuration 
in a wind tunnel at a fairly low Reynolds number 510,000. 

The following figures and plots, which represent both 
the theoretical and experimental results, are each shown 
with two sets of curves; one set represents the results 
with the winglets attached, while the other represents the 
results with no winglets. All configurations were tested at 
a gap and stagger both equal to one chord length. 

Figure 4-1 is a plot of the biplane lift coefficient 
versus angle of attack at a decalage angle of 0 degrees. 


The experimental results 

agree 

quite 

well 

with 

the 

theoretical 

results up 

to 

an angle 

of 

attack 

of 

approximately 

13 degrees. 

At this 

point 

the 

upper 

wing 

of 


the biplane begins to stall, while the lower wing remains 
unstalled until an angle of attack of approximately 18 
degrees is reached. Past 18 degrees, both wings are stalled 
and the lift coefficient drops off rapidly with increasing 
angle of attack. One of the primary reasons that biplanes 
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characteristically have a low maximum lift coefficient is 
because of the asymmetrical stall between the two wings. At 
a decalage angle of 0 degrees and a sjcagger of one chord 
length, the upper wing will operate at a higher lift 
coefficient than the lower wing. Theory predicts that at an 
angle of attack of approximately 12 degrees, the upper wing 

is operating at a C. of 0.924, while the lower wing is at 

L 

a C of only 0.556. This occurs because the lower wing 

1 

is emersed in the induced downwash of the upper wing, and 
hence, the lower wing operates at a less effective angle of 
attack. This effect does, however, bring on a smooth 
gradual stall rather than an abrupt stall. 

At the test Reynolds number of 510,000, the maximum 
lift coefficient for the no-winglet configuration was 
experimentally found to be 0.850. The same configuration 
tested with winglets showed a 3.6% increase in the maximum 
lift coefficient. It is expected that a much higher maximum 
lift coefficient would be reached in full scale flight 
tests. For example, an NACA 0012 two-dimensional 

symmetrical airfoil tested at a Reynolds number of 510,000 
has a maximum lift coefficient of aproximately 0.900 
(Ref. 19) where as the same airfoil tested at a Reynolds 
number of 6,000,000 yields a maximum lift coefficient of 


1 . 600 
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The biplane configuration tested with vinglets did show 
a slight increase in the lift-curve slope as would be 
expected. The theoretical lift-curve slope without winglets 
was computed to be 0.059 per degree, while the same 
configuration with winglets yielded a lift-curve slope of 
0.062 per degree; an increase of 5.1%. Experimentally, the 
lift-curve slope with w’inglets was 0.064 per degree, while 
the lift-curve slope without winglets was 0.061 per degree. 
These values were determined using algebraic linear 
regression. Also, as would be expected, the theoretical 
curves continue as a straight line since no flow separation 
is realized in potential flow theory. 

Figure 4-2 is a plot of the biplane lift coefficient 

versus angle of attack at a decalage angle of -5 degrees. 

At approximately 12 degrees angle of attack the experimental 

curve begins to diverge from the theoretical curve due to 

viscous flow separation. However, in this case the stall 

occurs much more abruptly due to both wings stalling at 

approximately the same time. Theory predicts that at an 

angle of attack of 12 degrees the upper wing is operating at 

a lift coefficient of 0.912 while the lower wing is at a 

lift coefficient of 0.953. Therefore, at this decalage 

angle a nearly symmetrical stall will occur between the two 

wings. For the no-winglet configuration, was 

max 

experimentally found to be 0.867 (which is slightly higher 

than for the no-winglet, zero decalage case). The 

max 
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biplane configuration (with a decalage angle of -5 degrees) 
tested with winglets showed a 4.1% increase in the maximum 
lift coefficient over the no-winglet case. 

The theoretical lift-curve slopes for the -5 degree 
decalage case, with and without winglets, were computed to 
be C.057 and 0.054 per degree respectively; a difference of 
5.3%. Experimentally, the lift-curve slopes were found to 
be 0.066 and 0.063 respectively. 


It can be reasoned that the change in decalage angle 
from 0 to -5 degrees has only a slight effect on the 
lift-curve slopes and the maximum lift coefficient. From 
figures 4-1 and 4-2 it is obvious that the effect of 
winglets is a slight increase in the maximum lift 
coefficient, as well as a slight increase in the lift-curve 
slopes . 


It is well known from previous wind tunnel testing that 
at low Reynolds numbers, profile drag coefficients can vary 
quite considerably with only relatively small variations in 
Reynolds numbers. Therefore, before the experimental data 
can be properly analyzed, it is desirable to know exactly 
what effects the low Reynolds numbers will have on profile 
drag coefficients. The profile drag coefficient, which 
consists of pressure (form) drag, skin friction, drag, , 
interference drag, and parasite drag, is commonly expressed 




+ * C D 

P 


a s : 


4.1 



Page 53 


where & is the incremental change in profile drag due to 
P 

lift. For the biplane -winglet configuration tested in this 
research, it was necessary to determine the relationship 
between the incremental profile drag coefficient and the 
lift coefficient at the test condition Reynolds number of 
510,000. This was done by studying the experimental results 
of Uenadovitch (Ref. 16), who in 1936 performed several 
experiments to determine the aerodynamic characteristics of 
two-dimensional biplane configurations utili 2 ing symmetrical 
airfoils. Figure 4-3 represents a plot of the incremental 
change in profile drag due to lift for a biplane 
configuration tested at a Reynolds number of 500,000. For 
this conf iguraion , the gap and stagger were both equal to 
one chord length. 


Figures 4-4 and 4-5 represent the experimental and 
theoretical drag polars for the biplane -winglet 
configuration tested at 0 and -5 degree decalage angles. 
The experimental points represent the data taken during 
wind-tunnel testing after the various wind-tunnel correction 
factors have been applied. The theoretical curve was 
plotted utilizing the following equation: 


+ L C. 


min 


4.2 


where is the minimum profile drag coefficient, which 

min 

for a symmetrical airfoil is also the profile drag 
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coefficient at zero lift. The values of the oiniirtum profile 

drag coefficients were determined from wind-tunnel tests. 

At a decalage angle of 0 degrees, the biplane configuration 

tested with winglets had a minimum profile drag coefficient 

of 0.017, while the same configuration without winglets 

yielded a minimum drag coefficient of 0.015. For the -5 

degree decalage case, the test results yielded minimum drag 

coefficients of 0.022 and 0.020, with and without winglets 

respectively. Thus, the presence of winglets adds an 

additional 0.002 to the minim urn profile drag coefficient at 

zero lift. The a term in equation (4.2) was determined 

P 

using figure 4-3, and the induced drag term, , was 

predicted theoretically using the vortex-lattice computer 
program. 

By referring to the drag polars in figures 4-4 and 4-5, 
the relative advantages of the winglets can be realized. At 
zero lift, the configuration with winglets produces slightly 
more drag (parasite drag caused by the presence of the 
winglets). However, as the lift coefficient begins to 
increase, the reduction in induced drag caused by the 
presence of the winglets begins to take affect. At lift 
coefficients greater than approximately 0.4, the 

configuration with winglets produces less total drag. This 
reduction in total drag becomes more significant as the lift 
coefficient increases further. For the -5 degree decalage 
case, the total drag is reduced by 6.5 % at a lift 
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coefficient of 0.7, while at a lift coefficient of 0.5 the 
drag is reduced by 3.3%. For the 0 decalage case, the 
advantages of adding winglets are not as significant as the 
-5 degree decalage case, however, above a lift coefficient 
of 0.4 there is still a reduction in the overall drag caused 
by the winglets. 

The experimental points plotted on the drag polars 
agree quite well with the theoretical curves, although the 
experimental points do show a consistently higher drag 
contribution. More importantly, the experimental points 
indicate approximately the same magnitude of drag reduction 
as do the theoretical curves for both configurations tested. 

It is important to realize that if these configurations 

were tested at a much higher Reynolds number, the induced 

drag reductions caused by the winglets would be greater. 

This is due to the significant effect of A at low 

P 

Reynolds numbers. For example, at an angle of attack of 8 
degrees, the biplane configuration with winglets produces a 
lift coefficient Of 0.519, which is a 6.6% increase over the 
lift coefficient produced by the same configuration with no 
winglets. This increase in lift coefficient causes the 
incremental profile drag coefficient to increase from 0.0144 
to 0.0166; an increase of 15.3%. At a Reynolds number of 
6,000,000 the profile drag coefficient will increase from 
0.0069 to 0.0072; an increase of only 4.3%. Therefore, 
full-scale advantages of winglets cannot be fully realized 
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at low Reynolds numbers. 

The biplane-wing let efficiency factors were computed by 

2 

first determining the various slopes of the C versus C 

L D 

curves. The efficiency factors were calculated based on the 
actual biplane aspect ratio, which is five. It is common in 
biplane theory to calculate the efficiency factors based on 
the biplane's equivalent monoplane aspect ratio, which for 
this configuration is 3.38. Using the equivalent monoplane 
aspect ratio would result in a 47% increase in the computed 
efficiency factors. 

Figures 4-6 and 4-7 represent the theoretical and 

2 

experimental plots of ' versus for the two decalage 

cases. The theoretical values were calculated by adding the 
minimum profile drag coefficient for a specific 

configuration to the theoretically predicted induced drag 
coefficients. The experimental points were determined by 
subtracting the incremental profile drag coefficients from 
the experimentally determined drag coefficients. The 
experimental drag coefficients then represent the minimum 
profile drag plus the induced drag. This was done in order 
to make logical comparisons between the theoretically 
predicted and experimentally determined biplane eficiency 
factors for the various configurations. 
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For the 0 decalage case, the theoretically predicted 
efficiency factors were found to be 0.737 and 0.683, with 
and without winglets respectively. Thus, the addition of 
winglets causes a theoretical increase in the efficiency 
factors of 7.9%. The efficiency factors were experimentally 
determined using algebraic linear regression techniques. 
The experimental values were determined to be 0.670 and 
0.588, with and without winglets respectively, for the 0 
decalage case. Thus, a 13.9% increase in efficiency was 
obtained experimentally by the addition of winglets. Also, 
it can be noted that the theoretical predictions are 
approximately 12% higher than the experimental values. 

For the -5 degree decalage case, the theoretical 
efficiency factors were found to to be 0.731 and 0.687, with 
and without winglets respectively. This is an increase of 
6.4%. Experimentally, the efficiency factors with and 
without winglets were calculated to be 0.663 and 0.562, a 
difference of 12.6%. Therefore, based on experimental 
results, the addition of winglets increases the efficiency 
factors by approximately 13%, which is a significant 
increase. The 0 degree decalage configuration was 
experimentally found to yield efficiency factors 
approximately 5% greater than the -5 degree decalage case, 
independent of whether winglets were attached or not. 
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Plots of the theoretically predicted induced drag 
coefficient as a function of the lift coefficient for the 
two decalage cases are shown in figures 4-8 and 4-9. It can 
be observed that as the lift coefficient increases, the 
reduction in induced drag afforded by the winglets becomes 
quite significant. At a lift coefficient of 0.6, the 
addition of winglets theoretically reduces the induced drag 
by 6.4%, while at a lift coefficient of 0.8, the induced 
drag is reduced by 8.3%. At a cruise lift coefficient of 
0.4 the induced drag is reduced by 6.2%. It is also 
significant to note that there appears to be little 
variation in the magnitude of induced drag reduction between 
the 0 and -5 degree decalage configuration. 



* u 



Figure 4-1 Lift Coefficient versus Angle of Attack 
for the 0 Degree Decalage Configuration 
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Figure 4-2 Lift Coefficient versus Angle of Attack 
for the -5 Degree Decalage Configuration 
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Figure 4-4 Biplane Drag Polar for the 0 Degree 
Decalage Configuration 
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Figure 4-5 Biplane Drag Polar for the -5 Degree 
Decalage Configuration 
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Figure 4-6 Lift Coefficient Squared versus Drag Coefficient 
for the 0 Degree Decalage Configuration 
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Figure 4-8 Theoretically Predicted Induced Drag ; t 
Coefficient versus Lift Coefficient for 
the 0 Degree Decalage Configuration 
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CHAPTER V 
COHCLl'S IOKS 

• t T • » ♦. : 

Based on results determined from the theoretical 
vortex-lattice computer program and the reduced experimental 
data collected from wind-tunnel tests conducted at a 

Reynolds number of 510,000, the following conclusions can 

can be determined concerning the aerodynamic characteristics 
of the biplane- winglet configuration. 

1. For the 0 degree decalage case, the addition of 

winglets increased the maximum lift coefficient by 
3.6%, increased the lift-curve slope by 5.1%, and 
.(based on experimental data) increased the overall 
efficiency factor by 13.4%. 

2. For the -5 degree decalage case, the addition of 

winglets increased the maximum lift coefficient by 
4.1%, increased the lift-curve slope by 5.3%, and 
increased the overall efficiency factor by 12.6% 

3. The most significant difference between the 0 and 
-5 degree decalage cases was that the 0 degree 
decalage case showed a 5% increase in the 
efficiency factor over the -5 degree decalage 
case. Also, the -5 degree decalage case showed a 
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slightly higher maximum lift coefficient# 


A. For both configurations, the addition of winglets 
resulted in approximately a 3.3% reduction in 
total drag at a lift coefficient of 0.5, and a 
6.5% reduction in total drag at a lift coefficient 
of 0.7. Below a lift coefficient of 0.A, the 


configurations 

tested 

with 

winglets produced 

si igh tly 

more total drag 

due 

to 

the additional 

parasite 

drag 

created 

by 

the 

winglets (the 


winglets produced an additional minimum profile 
drag increment of 0.002). 

5. The vortex-lattice computer program results showed 
a 6.2% reduction in induced drag at a lift 
coefficient of 0.A, and a 8.3% reduction at a lift 
coefficient of 0.8. Theoretical results did not 
indicate any significant differences in induced 
drag reductions (due to winglets) between the 0 
and -5 degree decalage cases. 


Overall , 

it can 

be 

determined that 

the 

addition 

of 

winglets 

to 

an already optimized 

biplane 

conf iguration 

can 

be beneficial 

with 

respect 

to 

increasing 

the maximum 

lift 

coef f icient , 

increasing 

the 

lift-curve slope, 

and 


increasing the overall efficiency of the lifting 
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system by decreasing the induced drag. 

Additional research is Suggested to further’ 
optimize the bi plane -winglet configuration. 
Possible extensions to the work already done are: 
varying the winglet airfoil and planform shape, 
and introducing large winglet cant and toe angles. 
These new configurations may further optimize the 
efficiency of the biplane-winglet configuration. 
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Append ix 


PROGRAM LISTING 



c 

f: 

C 

C: 

C 

C 

C 

c 

f: 

r: 

i.: 

C 

c 

c 

I* 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THIS PROGRAM WILL UTILIZE THE VORTEX LATTICE ME I HOD 
TO PREDICT THE AERODYNAMIC CHARACTERISTICS OF A 
B I PL ANE-W I NGLE T CONF I GURAT I ON . 

THE PRIMARY VARIABLES IN THIS PROGRAM WILL BE THE 
BIPLANE GAP, STAGGER, DECALAGE, SPAN, CHORD, ASPECT 
RATIO, WINGLET TOE ANGLE (INCIDENCE), WINGLET ASPECT 
RATIO, AND ANGLE OF ATTACK. 

FOR THIS PARTICULAR CONFIGURATION; 

VEL= FREESTREAM VELOCITY 
AR= ASPECT RATIO 
CH=CHORD 

GAP=GAP ( ONE CHORD LENGTH) 

ST AG=ST AGGER (IN PERCENT CHORD) 

DEC=DECALAGE ( DEGREES ) 

ALPHU= ANGLE OF ATTACK OF UPPER WING 
ALPHL= ANGLE OF ATTACK OF LOWER WING 
TOE=TC»E ANGLE OF WINGLET (POSITIVE- TOE INWARD) 

SW=WING PLANFORM AREA 

Nl= NO. OF HORSESHOE VORTICES LOCATED ON UPPER WING (HALF-SPAN) 
N2= NO. OF HORSESHOE VORTICES LOCATED ON BOTH UPPER 
AND LOWER WING (HALF-SPAN) 

N3= NO. OF HORSESHOE VORTICES LOCATED ON BOTH UPPER AND LOWER- 
WING PLUS THE WINGLET (HALF-SPAN) 


C 

C 

C 

C 

c 

c 

c: 


PN ( N ) , QN ( N ) , RN ( N ) = COORDINATES OF THE Nth HORSESHOE VORTEX 
PV(V) ,QV(V) ,RV(V)= COORDINATES OF THE Vth HORSESHOE VORTEX 

S= SEMI -WIDTH OF HORSESHOE VORTEX 


D I MENS I ON PN ( 1 00 ) , QN ( 1 00 ) , RN ( 1 00 ) , PV ( 1 00 ) , QV ( 1 00 ) , R V ( 1 00 ) 
DIMENSION FWS ( 1 00 , 1 00 ) , FWP ( 1 00 , 1 00 ) , F VS ( 1 00 , 1 00 ) , F VP ( 1 00 , 1 00 ) 
D I MENS I ON GAM ( 1 00 ) , WK ARE A ( 1 00 ) , A ( 1 00 , 1 00 ) , ^LPH ( 1 00 ) , B ( 1 00 ) 

D I MENS I ON PH ( 1 00 ) , WW ( 1 00 ) , WWL ( 1 00 ) , VW ( 1 00 ) , VWL ( 1 00 ) , TOE ( 1 00 ) 


« 
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DIMENSION SEC (100) 

INTEGER V. V1,V2,V3,V4,V5»V6, IA» I DOT 


PI=3. 141592 


VEL= 100.0 

**** INPUT BIPLANE-WINGLET GEOMETRY *#** 


WR I TE ( 6 » •* ) " ENTER 
READ (5. #)EW 
WR I TE ( 6 > # ) " ENTER 
READ (5. *)AR 
WRITE(6>» ■») ''ENTER 
READ (5,*) STAG 
WRITE( 6» * > ■'ENTER 
READ ( 5 » * ) ALPHU 
WR I TE ( 6 » * ) " ENTER 
READ (5,*) DEC 


WING TWIST" 

ASPECT RATIO" 
BIPLANE STAGGER" 
ALPHA (UPPER WING)" 
BIPLANE DECAL AGE" 


WRITE(6.*> "ENTER WINGLET MAXIMUM TOE ANGLE' 
RE AD ( 5 » # ) ALPH I 

BB=20. 0 

CH=BB/AR 

GAP=CH 

ALPHL=ALPHLI-DEC 

S=0. 25 

N1 =20 
N2=40 
N3=60 
N4=80 
N5=88 
N6=96 
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VI =20 
V2=40 
V3=60 
V4*80 
V5=8© 
V6=96 

KN=2 

K=8 


WR I TE < 8 , 200 )N6» BB» AR» ALRHU , DEC , STAG , CH , ALPH I 
200 FORMAT < // > 2X » 'TOTAL NUMER OF HORSESHOE VORTICES ■' , 1 4. / / . 2X , 'SPA 
GN = ' , F 4 . 1 » 25X » ' ASPECT RATIO » F&. 2» 15X » // » 2X » 'ALPHA UPPER ■ '» 
GF4.1.15X,' DECALAGE F4. 1 .//» 2X, 'STAGGER =',F4.1, 

G22X, 'CHORD ='»F4. 1 » //» 2X» 'MAXIMUM WINGLET TOE ANGLE =',F4.1> 


**** COMPUTE ARRAY OF COORDINATES FOR THE Nth HORSESHOE VORTEX **** 


DO 5 V=1 » N6 
DO 10 N=1,N6 

IF(N.LE.Nl) THEN 
PN(N)=— (CH/K) 

QN(N> = (2. *FLOAT <N)-1 . )*S 

RN<N)=0.0 

PH(N)=0.0 

ELSE IF(N.LE.N2> THEN 
PN < N ) =- ( 5. *CH/K ) 

QN(N)=( 2 . * ( FLOAT ( N ) -N l ) - 1 . ) *S 
RN<N)=O.C> 

PH ( N ) =0. O 

ELSE IFIN.LE.N3) THEN 
PN ( N ) =- ( STAG*CH ) - ( CH/K ) 

QN <N)=(2.»< N-N2 ) - 1 . )*S 

RN ( N ) =GAP 

PH(N)=0.0 
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ELSE I F ( N . LE . N4 ) THEN 

F'N ( N ) =- ( ST AG*CH ) - < 5#CH/K ) 

QN ( N ) = ( 2 . * < N-N3 ) — 1 . )*S 
RN ( N ) =GAP 
PH (IM) =0.0 

ELSE IF(N. LE. N5)THEN 

PN<N)=-<STAG*CH/<2.*K> >*<2. *<N-N4 )-l . )-(CH/K) 
QN(N)=BB/2. 


RN (N)=<2.*( N-N4 ) — 1 . )*(CH/16. ) 

PH(N)=PI/2. 

ELSE IF(N.LE.N6)THEN 

PN(N)=-(STAG*CH/<2.*K> >*<2. #(N-N5)-1 . )-<5. #CH/K) 

QN(N>=BB/2. 

RN(N)=<2.*(N-N5>-1. >*(CH/16. ) 

PH(N)=PI/2. 

END IF 

C **** COMPUTE ARRAY OF COORDINATES FOR THE Vth CONTROL POINT **** 

IF(V.LE.Vl) THEN 
PV(V)=-(3.*CH/K) 

Q V ( V ) = ( 2 . *FLOAT ( V ) — 1 . )*S 
RV(V)=0.0 
ALPH ( V ) =ALPHU 
ELSE IFIV.LE.V2) THEN 
PV ( V ) =— ( 7 . *CH/K ) 

QV(V) = <2.*(FL0AT(V)-V1)-1. )*S 
RV(V)=0.0 
ALPH ( V > =ALPHU 
ELSE IF (V.LE. V3) THEN 

PV < V ) — < STAG*CH ) - ( 3. *CH/K ) 

QV (V) = (2.*< V-V2 ) - 1 . ) *S 
RV(V)=GAP 
ALPH ( V ) =ALPHL 
ELSE IF(V.LE. V4)THEN 
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PV ( V ) =- < STAG*CH ) - ( 7 . *CH/K > 

Q V ( V ) = ( 2 . * ( V- V3 ) - 1 . ) *S 
RV(V)=GAP 
ALPH ( V ) =ALPHL 
ELSE IF (V.LE.V5) THEN 

PV<V)=-<STAG*CH/<2.*K> >*<2. *<V-V4)-1 . )-(3. *CH/K> 
QV< V)=BB/2. 

RV<V)=(2.*<V-V4)-1. >*<CH/16. ) 

TOE (V)=( GAP/2. -RV( V) >*<2. *< ALPHI/GAP) ) 

ELSE IF(V.LE. V6)THEN 

PV<V)=-<STAG#CH/<2.*K) )*<2. *< V-V5)-l . )-<7.*CH/K> 
QV ( V ) =BB/2. 

RV ( V ) = ( 2. # ( V-V5 ) - 1 . >*<CH/16. ) 

TOE ( V) =(GAP/2. — RV( V) >*<2.*< ALPHI/GAP) ) 

END IF 


X=PV( V)-PN(N) 

YS=QV(V)-QN(N) 

YP=QV(V)+QN(N) 

Z=RV(V)-RN(N) 

PHS*PH(N) 

PHP=-PH(N) 


*** COMPUTE DOWNWASH INFLUENCE COEFFICIENTS AT CONTROL POINTS ON WINGS*#* 

**** DOWNWASH COEFFICIENT- STARBOARD WING *** 

CALL FWSS ( FWS ,N,V,X,YS,Z,S, PHS , WS 1 , WS2 , WS3 , WS4 , WS5 , WS7 » WS9 ) 

**** DOWNWASH COEFFICIENT- PORT WING **** 

CALL FWPP ( FWP. N» V> X» YPt Z» S» PHP, WP1 , WP2» WPS, WP4 , WPS , WP7 , WP9 ) 
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c 

c 

c 

I" 

c 

c 

c 


*** COMPUTE SI DEWASH COEFICIENTS AT CONTROL POINTS ON MINGLE IS « * * 

*»»» SI DEWASH COEFFICIENT- STARBOARD MINGLE I ***# 

CALL FVSS ( F VS , N . V . X . YS . Z t. S . PHS , MS 1 , MS2 , WS3 » WS4 , MSS , MS7 , MS9 ) 

#*** SIDE WASH COEFFICIENT- PORT WINGLET **** 


CALL FVPP ( F VP , N , V . X , YP , Z , S , PHP , WP 1 , WP2 » WPS , MP4 , WPS , WP7 , WP9 ) 

C 

c 

C **# CALCULATE COEFFICIENT MATRIX FOR CONTROL POINTS LOCATED *»* 

C #** ON UPPER AND LOWER WINGS *#* 

C 

c 

IFtV.GT. V4)GO TO 100 

A ( V t N > =FWP ( V » N > +FWS ( V » N ) 

B(V>=4. #PI*VEL*ALPH(V)*<PI/180. ) 

GO TO 10 

0 

O *** CALCULATE COEFFICIENT MATRIX FOR CONTROL POINTS ON W1NGLETS 
C 

100 A(V,N)=FVP(V,N)+FVS(V,N) 

B< 81 >=4.*PI*VEL«TGE<81 >*<PI/180. ) 

B(82)=4. *PI*VEL#T0E<82)*(PI/180. ) 

B ( 83 ) =4 . *P I *VEL*TOE < 83 > # < P I / 1 80 . > 

B ( 84 ) =4 . #P I » VEL*T OE ( 84 ) * ( P I / 1 80. ) 
B(85)=4.*PI*VEL*T0E(35>*<PI/180. ) 

B(86)=4. *PI*VEL*T0E(86)*<PI/180. ) 

B ( 87 ) =4 . *P I *VEL*TOE ( 87 ) * ( P I / 1 80 . ) 
B<88)=4.*PI*VEL*T0E<88>#<PI/180. ) 

B(89)=B(8i ) 

B ( 90 ) =B ( 32 ) 

B<91 >=B<83) 

B(92)=B(84) 

B ( 93 ) =B ( 85 ) 
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B(94)=B(86> 

B ( 95 ) =B ( 87 ) 

B ( 96 ) =B ( 88 ) 

10 CONTINUE 
5 CONTINUE 

C 

C *** PROGRAM IS NOW READY TO SOLVE FOR IH£ GAMMA ARRAY FRUM THE *** 

C *** COEFFICIENT MATRIX AND THE BOIJNDRY CONDITION ARRAY *** 

C 

M=1 
N=N6 
I A= 1 00 
IDGT=6 

0 

0 *** CALL IMS'L SUBROUTINE TO SOLVE N SI MULTAN I OUS EQUATIONS #** 

C 

CALL LEQT 1F(A,M,N,IA,B.I DOT , WKAREA , I ER ) 

DO 15 V=1,N6 
GAM( V)=B( V) 

15 CONTINUE 


C #** CALCULATE LIFT COEFFICIENTS FOR UPPER AND LOWER WINGS *#* 

C 

CLUPP=0. 0 
CLL0W=0. O 

ST=2.*<BB*CH> 

SW=BB*CH 

DO 20 V= 1 , 20 
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CLUPP=CLUF'P+ ( 4 . / ( VEL*SW ) ) * ( ( GAM ( V ) +GAM ( V+20 ) ) #2 . 0*S ) 
SEC ( V ) = ( 2 . * ( GAM ( V ) +GAM ( V+20 ) > ) / ( CH*VEL ) 


1 1 0 CLLGW=CLLOW+ < 4 . / ( VEL*SW ) ) * ( ( GAM ( V+40 ) +GAM ( V+60 ) ) *2. 0#S ) 

SEC ( V+20 ) = ( 2 . * ( GAM ( V+40 ) +GAM ( V+60 ) ) > / < CH*VEL ) 

20 CONTINUE 

DO 21 V=8 1 * VS 

SEC ( V ) = ( 2. * ( GAM ( V ) +GAM ( V+S ) ) ) / ( CH*VEL ) 

21 CONTINUE 

CL- ( CLUPP+CLLOW )/ 2 . 

WRITE ( S» 18?) 

189 FORMAT (/.5X. 'CONTROL POINT' , 5X , 'SECTION LIFT COEFF' »3X, 'SECTION 
©LIFT COEFF', /, 27 X* 'UPPER WING' , 9X , 'LOWER WING',/) 

DO 16 V= 1 , 20 

WRITE (8, 190 )V, SEC( V) , SEC: (V+20) 

190 FORMAT ( 10X, 13, 12X,F9.5, 10X,F9.5) 

16 CONTINUE 

WRITE(8, 188) 

188 FORMAT < // , 5X, 'CONTROL POINT , 10X, 'WINGLET SECTION LIFT COEFF',/) 

DO 17 V=81 , V5 
WRITE (8, 1 87 ) V , SEC: ( V ) 

187 FORMAT ( 10X, 13, 22X , F9. 5) 

17 CONTINUE 

WR I TE < 8 , 202 ) CLUFP , CLLOW , CL 
202 F0RMAT(//,5X, 'LIFT COEFF (UPPER WING) =',F10.6, 

©/ , 5X, 'LIFT COEFF (LOWER WING) =',F10.6, 

@/,5X, 'TOTAL LIFT COEFF = ',F10.6) 

C 

C 

C **** CALCULATION OF INDUCED DRAG BY HUNK'S THEOREM OF INDUCED **** 
C ***# DRAG DETERMINATION UTILIZING THE TREFFTZ PLANE **** 
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TREF=-1 400. 0 


DO 50 V=1 i V6 
DO 55 N= 1 i N6 
X=TREF 

YS«QV(V)-QN(N) 
YP=QV ( V ) +QN ( N ) 
Z=R'V( V)-RN(N) 
F'HS=PH ( N ) 
PHP=-PH(N) 


CALL FWSS ( FWS ,N,V,X,YS,Z,S, PHS , WS 1 . WS2 » WS3 , WS4 , WS5 , WS7 , WS9 ) 


CALL FWPP(FWP,N, V, X, YP, Z.S.PHP.WPl , WP2, WP3, WP4, WPS, WP7.WP9) 
CALL FVSS ( FVS ,N,V,X,YS,Z,S, PHS , WS 1 , WS2 , WS3 , WS4 , WS5 , WS7 , WS*? > 
CALL FVPP ( F VP , N, V, X , YP» Z » S» PHP , WP 1 , WP2 , WPS , WP4 , WPS , WP7 , WP9 ) 
55 CONTINUE 
50 CONTINUE 


**** CALCULATION OF INDUCED VELOCITYS IN THE TREFFTZ PLANE *»#* 


R0W=0. 002378 
DO 60 V=1 » 96 


WW(V)=0.0 
WWL. ( V ) =0.0 
VW(V)=0.0 
VWL ( V ) =0 . O 

DO 65 N= 1 » 80 

I F ( V . LE . 80 ) THEN 

WW < V ) = ( GAM ( N ) / ( 4 . #P I ) )*<FWP(V,N)+FWS<V,N) )+WW(V) 
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ELSE I F ( V . LE . 96 ) THEN 

VW( V)=<GAM<N)/ (4. #PI ) )*(FVP(V,N)+FVS< V,N) >+VW<V) 

END IF 

65 CONTINUE 

DO 70 N=81 » 96 

I F < V . LE . 80 ) THEN 

WWL<V>=(GAM(N>/<4. *FI > > * ( FWP < V » N ) +FWS ( V , N ) )+WWL.(V> 

ELSE I F ( V . LE . 96 ) THEN 

VWL ( V ) = (GAM ( N ) / ( 4 . #PI ) > * ( FVP< V, N) +FVS ( V, N ) )+VWL< V) 

END IF 

7u CONTINUE 

60 CONTINUE 

C 

c 

C *** CALCULATION OF INDUCED DRAG CREATED BY UPPER AND LOWER WINGS *** 

C 

C 

DWW=0. 0 
DWWL=0. 0 
DVW=0.0 
DVWL=0. O 


DO 75 V= 1 , SO 
N=V 

DWW=R0W*2. *S* WW ( V ) *GAM ( N ) +DWW 
DWWL=R0W*2. *S*WWL ( V ) *GAM ( N > +DWWL 
75 CONTINUE 

C 

C 

C *#* CALCULATION OF INDUCED DRAG CREATED BY WINGLETS *** 
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DO 80 V=81,96 
N=V 

DVW=-R0W*2. *S*VW ( V ) *GAM ( N ) +DVW 
DVWL=-R0W*2 . *S*VWL ( V ) *GAM ( N ) +DVWL 
80 CONTINUE 

C 

c 

C *#* CALCULATION OF DRAG COEFFICIENTS FOR WINGS AND WINGLETS 

C 

C 

0=0. 5#R0W* ( VEL#*2. ) 

CDWW=DWW/(Q*ST) 

CDWWL=DWWL/ ( Q#ST ) 

CDVW=DVW/(Q*ST> 

CDVWL=DVWL/<0*ST) 

CDW=CDWW+CDWWL 

CDWL=CDVW+CDVWL 

CD=CDW+CDWL 

E= ( CL**2 . O > / < P I *AR*CD > 

CLALPH=CL/ ( ALFHU- ( DEC/ 1 . 667 ) ) 

RLD=CL/CD 

WR I TE ( 8 , 203 ) CDWW , CDWWL , CDVW , CDVWL , CDW , C DWL , CD» E. CLALPH > RLD 
203 FORMAT ( / / , 5X , ' DRAG COEFF ( W I NG-W I NG ) = , F 1 6 . 6 , 

0/ > 5X * •'DRAG COEFF ( W I NG-W I NGLET ) = •' » F 1 3 . 6 > 

© / , 5 X , ' DRAG COEFF < W I NGLET -W I NO > = ' , F 1 3 . 6 , 

© / , 5 X » ' DRAG COEFF ( W I NGLET-W I NGLET > = ' , F 1 0 . 6 , 

0//.5X, 'TOTAL WING DRAG COEFF= ' , F13. 6, 

0/.5X, 'TOTAL WINGLET DRAG COEFF=' , F10. 6, 

0/,5X, TOTAL BIPLANE DRAG COEFF=' , F10. 6, 

0 / / > 5 X , • OSWALDS EFF I C I ENC Y FACTOR= ' , F 1 0 . 6 , 

0 / i 5X , ' 3-D L I FT SLOPE CURVE ( /DEG ) = , F9 . 6 , 

0/.5X, LIFT/DRAG RAT 10 (W/ CDmin>= ,F9.6) 


♦ 
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/ > 


-5 


***** THE following section of this program will compute the ***** 
* * 
***** AERODYNAMIC CHACTERISTICS UF THE BIPLANE CONFIGURATION ***** 
* * 
******** ***tt»*»»«tt«* WITH NO WINGLETS **************************** 


DO 6 V=1 » N4 
DO 11 N=1 > N4 


X =PV ( V ) -PN ( N > 
YS=QV(V)-QNTN> 
YP=QV(V)+QN(N) 
Z=RV( V)-RN(N) 
PHS=PH(N> 
F'HP=-PH ( N ) 


*** COMPUTE DOWNWASH INFLUENCE COEFFICIENTS AT CONTROL POINTS ON WINGS*** 


**** DOWNWASH COEFFICIENT- STARBOARD WING *** 


CALL FWSS <FWS» N»V*X»YS*Z»S» PHS , WS 1 » WS2 , WS3 » WS4 » WS5 , WS7 » WS9 > 
**** DOWNWASH COEFFICIENT- PORT WING **** 

CALL F WPP ( FWP » N » V > X » YP » Z > S » PHP , WP 1 , WP2 » WPS , WP4 , WP5 » WP7 » WP9 ) 


*** CALCULATE COEFFICIENT MATRIX FOR CONTROL POINTS LOCATED *** 

*** ON UPPER AND LOWER WINGS *** 



A< V,N)=FWP(V,N)+FWS(V,N) 

B<V)=4. *PI*VEL#ALPH<V)*<PI/180. ) 

11 CONTINUE 
6 CONTINUE 

PROGRAM IS NOW READY TO SOLVE FOR THE GAMMA ARRAY FROM I HE 
**» COEFFICIENT MATRIX AND THE BOUNDRY CONDITION ARRAY *** 


M=1 
N=N4 
I A=1 00 
IDGT=6 

**» CALL I MSL SUBROUTINE TO SOLVE N SI MULTAN I OUS EQUATIONS *** 

CALL L.EQT 1 F ( A . M > N » I A » B > I DOT. WK AREA » I ER ) 

DO 26 V=1,N4 
GAM( V)=B( V) 

6 CONT I NUE 


#*# CALCULATE LIFT COEFFICIENTS FOR UPPER AND LOWER WINGS *** 

CLUF'P=0. 0 
CLL0W=0. O 

ST=2. * < BB#CH ) 

SW=BB*CH 

DO 22 V= 1 » 20 

CLUPP=CLUPP+ ( 4 . / ( VEL*SW ) ) * ( ( GAM ( V ) +GAM ( V+20 ) )*2. 0*S ) 

SEC ( V ) = ( 2 . * < GAM ( V ) +GAM ( V+20 ) ) ) / ( CH*VEL ) 

CLLOW=CLLOW+ ( 4 . / ( VEL*SW ) ) * ( ( GAM ( V+40 ) +GAM ( V+60 ) ) *2 . 0*S ) 
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t'j «:n 


♦ 


SEC ( V+20 ) = ( 2. * ( GAM < V+40 ) +GAM ( V+60 ) ) ) / ( CH*VEL ) 

22 CONTINUE 

CL= < CLUPP+CLLOW ) /2. 

WRITE! 8, 149 ) 

149 FORMAT!///, TOTAL NUMBER OF HORSESHOE VORTICES =',I4,5X, 

@ •' NO WINGLETS! ", //> 

WRITE (8. 150) 

150 FORMAT !/,5X, CONTROL POINT 5X, SECTION LIFT COEFF', 3X, SECTION 
©LIFT COEFF' » / , 27X . •'UPPER WING' ,9X, LOWER WING',/) 

DO 27 V-1.20 

WR I TE ! 3 » 151 )V, SEC! V), SEC! V+20) 


1 FORMAT! 10X, 13. 12X.F9.5, 1QX.F9.5) 
7 CONTINUE 


WRITE! 8, 152)CLUPP, CLLOW, CL 
152 FORMAT ! // , 5X , •'LIFT COEFF ! UPPER WING) =',F10.6, 

©/,5X, LIFT COEFF ILOWER WING) =',F10.6, 

©/ * 5X » •'TOTAL LIFT COEFF = '.F10.6) 

C 

C 

C **** CALCULATION OF INDUCED DRAG BY HUNK'S THEOREM OF INDUCED ###* 
C **»* DRAG DETERMINATION UTILIZING THE TREFFTZ PLANE **** 

C 

C 

TREF=- 1400.0 

DO 51 V=1,V4 
DO 56 N=1,N4 

X=TREF 

YS=QV(V)-QN!N) 

YP=QV ( V ) +QN ! N ) 
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Z=RV( V)-RN(N) 

PHS=PH ( N ) 

PHP=-PH(N) 

CALL FUSS ( FWS , N , V , X , YS , Z > S , PHS , WS 1 , WS2 , W.S3 , WS4 , WS5 , WS7 , WS9 ) 
CALL FWPP ( FWP .N>V>X» YP» Z» S> PHP , WP 1 , WP2 , WPS , WP4 , WPS , WP7 , WP9 ) 
CALL FVSS < F VS ,N,V,X>YS,Z,S, PHS , WS 1 , WS2 , WS 3 > WS4 , WS5 , WS7 , WS9 ) 
CALL FVPP ( FVP ,N»V,X,YP,Z,S, PHP , WP 1 , WP2 , WPS , WP4 , WP5 , WP7 , WP9 » 

56 CONTINUE 
51 CONTINUE 

C 

C 

C **#* CALCULATION OF INDUCED VELOCITYS IN THE TREFFTZ PLANE #*** 
C 

c 

R0W=0. 002378 

DO 61 V= 1 « 80 

WW<V)=0.0 
WWL ( V ) =0.0 
VW< V)=0.0 
VWL ( V ) =0 . O 

DO 66 N= 1 > 80 

I F ( V . LE . 80 ) THEN 

WW ( V ) = ( GAM <N)/(4.»PI > >*<FWP<V»N)+FWS<V,N> >+WW<V> 

ELSE I F ( V . LE . 96 ) THEN 

VW ( V > = ( GAM ( N ) / ( 4 . *P I ) ) * ( FVP ( V » N ) +F VS < V » N ) )+VW< V) 

END IF 

66 CONTINUE 

61 CONTINUE 

C 

o 

e 
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•» # » 


c ««* calculation of induced drag created by upper and lower minus 

o 


c 

DWW=0. 0 
DWWL=0. O 
DVW=0. 0 
DVWL=0.0 


DO 76 V=i,80 
N=V 

DWW=R0W*2 . *S*WW ( V ) *GAM ( N ) +DWW 
DWWL=R0W*2 . *S*WWL ( V ) *GAM ( N ) +DWWL 
76 CONTINUE 

C 

c 

C *** CALCULATION OF DRAG COEFFICIENTS FOR WING **# 

C 

i": 

G=0.5*R0W*<VEL**2. ) 

COWW=DWW/(Q*ST) 

CDWWL=DWWL/ ( Q*ST ) 

CDVW=DVW/(Q*ST) 

CDVWL=DVWL/ ( Q*ST ) 

CDW=CDWW+CDWWL 

CDWL=CDVW+CDVWl_ 

CD=CDW+CDWL 

E== ( CL**2. O) / ( PI *AR*CD) 

CLALPH=CL/ ( ALPHU- ( DEC/ 1 . 667 ) ) 

RLD=CL/CD 

WR I TE < 8 , 1 54 > CDMM , CDWWL > CDVW » CDVWL » CDW » CDWL , CD , E , CL ALPH , RLD 
1 54 FORMAT ( / / , 5X , DRAG COEFF ( W I NG-W I NG > = •' , F 1 6 . 6 , 

@/,5X, DRAG COEFF( WING-WINGLET ) = •' » F13. 6 * 
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®/,5X, DRAG COEFF ( W I NGLET-W I NO ) = ' > F 1 3 . 6 » 
@/.5X, "DRAG COEFF ( W JNGLET-WINGLET ) = " , FI 0. 6, 
@//.5X, "TOTAL WING DRAG COEFF= ' , F 13. 6 , 
@/,5X. "TOTAL WINGLET DRAG COEFF= " , F10. 6, 
@/,5X, "TOTAL BIPLANE DRAG COEFF=" , F 10. 6, 

@/ / > 5X » OSWALDS EFF I C I ENCY FACTOR= " , F 1 0. 6 , 
@/.5X, '3-D L I F T SLOPE CURVE < / DEG )=" ,F9. 6, 

©/ 1 5X > 'LIFT /DRAG RATIO! W/ CDmir.) = ",F9.6> 


1000 STOP 
END 


C ************* SUBROUTINE FWS *************** 
C * * 
C **» DOWNWASH COEFFICIENT- STARBOARD WING *** 


SUBROUT I NE FWSS ( FWS , N , V , X , YS » Z . S , PHS , WS 1 , WS2 , WS3 , WS4 » WS5 , WS7 , WS9 ) 
D I MENS I ON FWS < 1 00 .100) 

INTEGER V 

WS 1 = ( - X *C:QS < F’HS ) )/( ( X**2. ) + < ( Z*COS(PHS ) -YS*SIN ( PHS ) )**2. ) ) 

WS2= ( < YS+S*COS ( PHS ) ) *COS ( PHS ) ) + ( ( Z +S*S I N < PHS ) ) *S I N ( PHS ) ) 

WS3= ( ( X**2. )+( YS+S*COS(PHS) >**2. +< Z+S*SIN(PHS) )**2. >**0.5 
WS4= ( ( YS-S*COS < PHS ) )*COS<PHS> >-M < Z-S#SIN<PHS> )*SIN(PHS) ) 

WS5=< ( X**2. )+< YS-S*COS(PHS) )**2. +( Z-S*SIN(PHS) )**2. )**0.5 
WS6= ( YS-S*COS ( PHS ) ) / ( ( YS-S*COS ( PHS ) )**2. +( Z-S*SIN(PHS) )**2. ) 
WS7= 1 . - ( X/WS5) 


WS8= ( YS+S*COS ( PHS ) )/( ( YS+S*COS(PHS) ) **2. +( Z+S*SIN(PHS> )**2. ) 
WS9=1 . - ( X/WS3) 

FWS ( V , N ) =WS 1 * ( WS2/WS3-WS4/WS5 ) - ( WS6*WS7 ) + ( WS8*WS9 ) 

RETURN 

END 




V 
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SUE«ROUT INE FWP #*# ## ***** # 

* # 

*** DOWNWASH COEFFICIENT- PORT WING *** 

SUBROUT I NF. FWPP ( FWP , N , V , X , YP , Z , S , F'HP , WP 1 , WP2 , WPS , MP4 > WPS > WP7 , MP9 ) 
D I MENS I ON F WP ( 1 00 > 1 00 ) 

INTEGER V 

WP 1 = < - X *COS ( PHP ) )/( <X**2. > + < (Z*COS(PHP)-YP*SIN(PHP) )**2. ) ) 

WP2= ( ( YP+S*COS ( PHP ) ) #COS < PHP ) ) + ( ( Z +S*S I N < PHP ) ) * S I N < PHP ) ) 

WP3= ( (X**2. ) + ( YP+S*COS(PHP) )##2. +< Z + 8« i- IN (F'HP) )##2. >**0.5 
WP4= ( ( YP— S*COS ( PHP ) ) *CQS ( PHP > ) + ( ( Z-S*S I N ( PHP ) ) *S I N ( PHF ) ) 

WP5=< <X**2. )+< YP-S*COS(PHP) >**2. +< Z-S»SIN<PHP) >**2. >**0.5 
WP6= ( YP-S»CGS ( PHP > )/( ( YP-S*COS(PHP) ) **2. + ( Z-S*S IN ( PHP ) >**2. ) 
WP7=1 . -( X/WP5) 

WP8= ( YP+S*COS ( PHP ) )/( ( YP+S*COS(PHP> > **2. + ( Z+S*S1N ( PHP ) )**2. ) 
WP9= l . — ( X/WP3 ) 

FWP ( V, N ) =WP1 * ( WP2/WP3-WP4/WP5) - ( WP6*WP7 ) + ( WP8*WP9 > 

RETURN 

END 


**«*****#♦#«* SUBROUTINE FVS **##«***#*##**#### 

* * 

*** SIDEWASH COEFFICIENT- STARBOARD WINGLET *** 

SUBROUT I NE FVSS ( FVS >N,V,X,YS,Z,S, PHS , WS 1 , WS2 , WS3 . WS4 , WS5 , WS7 , WS 
D I MENS I ON FVS ( 1 00 , 1 00 > 

INTEGER V 

VS1 = < X*SIN(PHS) )/( (X**2. >+<Z*CQS(PHS>-YS*SIN(PHS> >**2. > 

VS2=WS2 

VS3=WS3 

VS4=WS4 

VS5-WS5 

VS6= ( Z-S*SIN(F'HS ) )/( ( YS-S*COS < PHS ) ) **2. + < Z-S*SIN ( PHS ) >**2. ) 
VS7=WS7 

VS8=(Z+S*SIN(PHS> )/ ( (YS+S*COS(PHS> >**2.+(Z+S*SIN(PHS> >**2. > 
VS9=WS9 


9) 


FVS ( V , N ) =VS 1 * ( VS2/ VS3-VS4/ VS5 ) + ( VS6*VS7 ) - ( VS8*VS9 ) 
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RETURN 

END 


C ####***##** SUBROUTINE FVP #«■*#*#*#**####* 

C * ' # 

C *** SI DEWASH COEFFICIENT- PORT WINGLET **# 

C 

SUBROUT I NE FVPP ( FVP , N , V , X , YP , Z , S , PHP , WP 1 , WP2 , WPS , WP4 » WPS , WP7 , WP9 ) 
D I MENS I ON FVP ( 1 00 , 1 00 ) 

INTEGER V 

VP 1 = ( X *S I N ( PHP ) ) / ( (X#*2. >+<Z#COS(PHP)-YP*SIN<PHP) )*#2. ) 

VP2=WP2 

VP3=WP3 

VP4=WP4 

VP5=WP5 

VP6= ( Z -S#S I N ( F'HP ) )/( ( YP-S*COS(PHP> )**2. +( Z-S*SIN(PHP) )**2. ) 
VP7=WP7 


VP8=<Z+S*SIN<PHP> )/( <YP+S*COS<PHP> )»*2. + ( Z+S*SIN < PHP ) >**2. ) 
VP9=WP9 

F VP < V , N ) =VP 1 * ( VP2/VP3-VP4 / VPS ) + ( VP6»VP7 ) - ( VPS# VP9 ) 

RETURN 

END 




4 
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